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In this paper we explore the effect of decaying dark matter (DDM) on large-scale structure and 
possible constraints from galaxy imaging surveys. DDM models have been studied, in part, as a 
way to address apparent discrepancies between the predictions of standard cold dark matter models 
and observations of galactic structure. Our study is aimed at developing independent constraints 
on these models. In such models, DDM decays into a less massive, stable dark matter (SDM) 
particle and a significantly lighter particle. The small mass splitting between the parent DDM and 
the daughter SDM provides the SDM with a recoil or "kick" velocity vu, inducing a free-streaming 
suppression of matter fluctuations. This suppression may be probed via weak lensing power spectra 
measured by a number of forthcoming imaging surveys that aim primarily to constrain dark energy. 
Using scales on which linear perturbation theory alone is valid (multipoles £ < 300), surveys like 
Euclid or LSST can be sensitive to Vk > 90 km/s for lifetimes r ~ 1 — 5 Gyr. To estimate more 
aggressive constraints, we model nonlinear corrections to lensing power using a simple halo evolution 
model that is in good agreement with numerical simulations. In our most ambitious forecasts, using 
multipoles £ < 3000, we find that imaging surveys can be sensitive to Vk ~ 10 km/s for lifetimes 
t < 10 Gyr. Lensing will provide a particularly interesting complement to existing constraints in 
that they will probe the long lifetime regime (r S> Hq 1 ) far better than contemporary techniques. 
A caveat to these ambitious forecasts is that the evolution of perturbations on nonlinear scales will 
need to be well calibrated by numerical simulations before they can be realized. This work motivates 
the pursuit of such a numerical simulation campaign to constrain dark matter with cosmological 
weak lensing. 

PACS numbers: 95.35.+d,98.80.-k,98.62.Gq, 
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I. INTRODUCTION 

Many and various astronomical observations indicate 
that ~ 5/6 of the mass density of the Universe is non- 
baryonic dark matter (reviews include Refs. [l|-|3| ) . The 
simplest model of so-called cold dark matter (CDM) can 
be successfully applied to interpret an enormous amount 
of observational data, particularly those characterizing 
the large-scale (> a few Mpc) structure of the Universe 
and the gross properties of galaxies. In particular, the 
CDM model is consistent with the cosmic microwave 
background (CMB) anisotropy spectrum measured by 
the Wilkinson Microwave Anisotropy Probe (WMAP) 
and observations of the large-scale (k < 0.1 h/Mpc) 
galaxy clustering spectrum measured by the Sloan Digital 
Sky Survey (SDSS) [J]. On smaller scales, the situation is 
murkier. Several observations indicate possible discrep- 
ancies between CDM theory and observations on smaller 
scales. Among these are the well-known missing satellites 
problem and the steep rotation curves of low-surface 
brightness galaxies [7- 9] . Exploring small-scale structure 
is challenging both observationally and theoretically. On 
the theoretical side, it is necessary to model highly non- 
linear phenomena to predict the properties of galaxies 
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and the dark matter halos in which they reside. Never- 
theless, these potential shortcomings of CDM may point 
toward novel properties of dark matter and many alterna- 
tives to CDM have been considered, including warm dark 
matter (WDM) [IMl and self-interacting dark matter 
(SIDM) [H-[l7|]. 

As one alternative to CDM, unstable dark matter has 
been considered in a number of recent studies (l5l 
[2f| . In such models, a dark matter particle of mass 
M decays into a less massive daughter particle of mass 
rn = (1 — /) M and a significantly lighter, relativistic 
particle, with a lifetime on the order of the age of the 
Universe. In decaying dark matter models, cosmolog- 
ical structure growth is altered in a time- and scale- 
dependent manner [2(| We explored models with 
/ ~ 1 in a previous paper (26|. If / <§; 1, the sta- 
ble dark matter (SDM) daughter particle will receive a 
non-relativistic kick velocity, Vk — f c. This kick can be 
sufficient to alter small-scale structure growth, modify- 
ing the predicted structures and abundances of the dark 
matter halos that host galaxies. Previous studies have 
explored just this possibility and placed limits on these 
scenarios by demanding that the alterations to structure 
growth not be so severe as to destroy the successes that 
CDM theory enjoys on large scales. An upper bound 
on the lifetime of the decaying dark matter (DDM) par- 
ticle of r > 30 — 40 Gyr for Vk ^ 100 km/s can be de- 
rived from the observed galaxy-cluster mass function and 
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galaxy mass-concentration relation [ID, [l7[. Lifetimes 
r < 30 Gyr with 20 km/s < < 200 km/s may be in 
tension with the observed properties of the Milky Way 
satellite galaxy population, but uncertainties associated 
with the details of nonlinear structure growth in these 
models as well as the star formation histories of Local 
Group galaxies are significant [22j]. Significantly tighter 
constraints can be obtained if the daughter particles are 
standard model particles [23] ]. 

With a number of forthcoming large-scale galaxy sur- 
veys being undertaken this decade, such as the Dark 
Energy Survey (DES) 1 , the Panoramic Survey Telescope 
and Rapid Response System (Pan STARRS) 2 , the Large 
Synoptic Survey Telescope (LSST) 3 |3, Euclid 4 [28|, 
and WFIRST 5 , it will be possible to investigate possi- 
ble subtle imprints of decaying dark matter on the rel- 
atively large-scale structure of the Universe. In DDM 
models, the kick velocity at decay imparts upon the sta- 
ble daughter particles the ability to smooth gravitational 
potential perturbations on scales smaller than the classic 
free-streaming scale. This behavior is similar to the cos- 
mological influence of massive neutrinos or WDM, and 
numerous studies have shown that these features may be 
detectable in large-scale structure through galaxy cluster- 
ing [29| , Lyman-q forest [30T - [32j , and cosmological weak 
lensing |33l - [35| . 

In this paper, we explore the effect of DDM model 
on large-scale structure and study the possible indepen- 
dent constraints on these models from forthcoming weak 
lensing surveys. There are distinct advantages to this ap- 
proach. First, DDM models are explored largely in order 
to mitigate the possible small-scale problems of CDM, 
so it is necessary to explore independent predictions and 
probes in order to test such models. Second, the effects of 
DDM on large-scale structure, and on weak lensing power 
spectra in particular, are simpler to model in that they 
do not require detailed knowledge of galaxy formation in 
the highly nonlinear regime including star formation his- 
tories, star formation and active galaxy feedback, scale- 
dependent galaxy clustering bias and numerous other 
complicated phenomena that are known to be important 
on small scales. Third, many surveys already have the 
goal of measuring cosmological weak lensing as a probe 
of dark energy 36] , so this test can be performed largely 
with the observational infrastructure used to study dark 
energy at no additional cost. Lyman-a forest spectra pro- 
vide an additional, promising method to constrain DDM, 
but also introduce complications associated with model- 
ing the clustering of neutral hydrogen along lines of sight 
to high redshift. We will present Lyman-a constraints 
on DDM, and forecasts for future measurements of the 
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Lyman-a forest, in a forthcoming follow-up paper. 

In the following sections, we will show that the im- 
print of DDM on lensing power spectra is sufficiently dis- 
tinct from other cosmological parameters, such as neu- 
trino mass, that one can disentangle degeneracies among 
them. Indeed, forthcoming data will be able to disentan- 
gle the two parameters of DDM models, lifetime and kick 
velocity (or equivalently /). Finally, large-scale lensing 
surveys will provide, at minimum, competitive and in- 
dependent constraints on DDM models exploiting only 
scales on which linear cosmological perturbation theory 
is appropriate. These probes will have the particular 
advantage of probing significantly larger DDM lifetimes 
than small-scale structure studies and will not require 
detailed modeling of galaxy formation. Moreover, a rela- 
tively modest numerical simulation program may enable 
one to use mildly nonlinear scales fc < 1 h/Mpc to obtain 
constraints on DDM that exceed current constraints and 
are robust to uncertainties associated with star formation 
and feedback in small galaxies. 

The outline of this paper is as follows. In § HH we de- 
scribe cosmological weak lensing observables. In § IIII1 
we discuss the perturbation evolution of both the par- 
ent and daughter particles in DDM models. In § IIV1 
we study the effect of DDM on dark matter halo den- 
sity profiles and explore possible influences of DDM on 
nonlinear structure. In § |Vl we describe the methods we 
use to compute constraints on DDM model parameters. 
We present our results in § IVI[ beginning with the gen- 
eral effects of DDM on large-scale structure, including a 
detailed discussion of the free-streaming of the daughter 
SDM particles. We also give our forecasts for constraints 
on DDM from large-scale weak lensing surveys and com- 
pare possible future bounds with existing limits. In § IVIII 
we summarize our work. 



II. WEAK GRAVITATIONAL LENSING 
OBSERVABLES 

Weak lensing as a cosmological probe has been dis- 
cussed at length in numerous papers (a recent review is 
Ref. We give a brief description of our methods be- 

low, which are based on the conventions and notation in 
Ref. [38| . The most robust forecasts derive from consid- 
erations of possible weak lensing measurements restricted 
to scales where linear perturbative evolution of the metric 
potentials remains useful. However, we attempt to esti- 
mate possible improvements to the constraining power of 
weak lensing observables, provided that mildly nonlinear 
evolution can be modeled robustly. 

We consider the set of observables that may be avail- 
able from large-scale galaxy imaging surveys to be the 
auto- and cross-spectra of lensing convergence from sets 
of galaxies in iVroM redshift bins. The ./Vtom (-^tom + 
l)/2 distinct convergence spectra are 

™ = / dz Zt)DitJ) Pm{k = e/DA > z) > (1) 
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where i and j label the redshift bins of the source galax- 
ies. We take Atom = 5 an d evenly space bins in red- 
shift from a minimum redshift of z = to a maximum 
redshift of z — 3. Increasing the number of bins beyond 
Atom = 5 adds only negligibly to the constraining power 
of lensing data, in accord with an analogous statement 
for dark energy constraints [39(. 

In Eq. ([1]), H(z) is the Hubble expansion rate, Da(z) 
is the comoving angular diameter distance, and P m (k, z) 
is the matter power spectrum at wavenumber k and red- 
shift z. In the following section, we describe our use of the 
publicly- available CMBFAST code to calculate P m (k) from 
cosmological perturbation evolution. In this case it will 
be more natural to work in the synchronous gauge, and 
transforming between different gauge systems can be ac- 
complished straightforwardly by following, for example, 
the methods described in Ref. |40l |. 

The Wi are the lensing weight functions for source 
galaxies in redshift bin i. In practice, the galaxies will be 
binned by photometric redshift, so that the bins will have 
non-trivial overlap in true redshift (see Refs. [HI, |4l| for 
detailed discussions) . Defining the true redshift distribu- 
tion of source galaxies in the zth photometric redshift bin 
as drii/dz, the weights are 



Wi(z) 



n M H 2 (l + z)D A (z) J dz 



, Da{z, z') dn t 
D A (z') Iz 7 



(2) 



where Da(z, z') is the angular diameter distance between 
redshift z and z' and Hq is the present Hubble rate. 

We model the uncertainty induced by utilizing photo- 
metric galaxy redshifts with the probability function of 
assigning an individual source galaxy photometric red- 
shift z p given a true redshift z, P{z p \z). The true redshift 
distribution of sources in the ith photometric redshift bin 



drii(z) 



dz r 



dn(z) 



P(z p \z) 



dz J J)ow) dz 

Here we take the true redshift distribution to be 

dn(z) _ 4z 2 



dz 



V27 



: exp[-(z/z ) 



(3) 



(4) 



with Zq — 0.92, so that the median survey redshift to 
Zmed = 1) an d n as the total density of source galaxies 
per unit solid angle [42Tl44j |. We assume that uncertain 
photometric redshifts can be approximated by taking 



P(z p \z) 



1 



V2trjI 



exp 



(z p - zf 



2(7? 



(5) 



where a z (z) = 0.05(1 + z) [39|. Complexity in photomet- 
ric redshift distributions is an issue that must be over- 
come to bring weak lensing constraints on cosmology to 
fruition (e.g., Ref. jHEI). 

Observed convergence power spectra P^(£), contain 
both signal and shot noise, 
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where (j 2 ) is the noise from intrinsic ellipticities of source 
galaxies, and rii is the surface density of galaxies in the 
ith tomographic bin. We follow recent convention and 
set yj (7 2 ) =0.2, subsuming additional errors on galaxy 
shape measurements into an effective mean number den- 
sity of galaxies, n. Assessments of intrinsic shape noise 
per galaxy may be found in, for example [13, HH, E3] . 
Assuming Gaussianity of the lensing field, the covariance 
between observables P l J and P. M is 



Cai 



pik pjl _j_ pit pjk 



(7) 



where the i and j map to the observable index A, and 
k and 1 map to B such that Cab is a square covariance 
matrix with ./Vtom(Atom + l)/2 rows and columns. We 
assume Gaussianity throughout this work and even in 
our most aggressive forecases we consider only multipoles 
£ < 3000, at which point the Gaussian assumption and 
several weak lensing approximations break down [48|-[52j . 



III. DECAYING DARK MATTER MODELS 

A. Evolution of the Average Properties of 
Unstable Dark Matter 

We consider dark matter decays into another species 
of stable dark matter with a small mass splitting, 
DDM — > SDM + L, where L denotes a "massless" daugh- 
ter particle, SDM is the stable dark matter with mass m, 
and DDM is the decaying dark matter with mass M . The 
mass loss fraction f of DDM is directly related to the kick 
velocity deposited to the SDM particle by / ~ /c from 
energy-momentum conservation. The following relations 
are valid in the rest frame of DDM particles with the kick 
velocity of SDM being the velocity relative to the DDM 
rest frame. 

For a general decay, neglecting Pauli-blocking factors 
and inverse decays, the rate of change in the DDM dis- 
tribution function is 



j ddm (qd dm) — — 



aMT 
Eddm 



Jddm {qddm), (8) 



where / denotes the partial derivative of the distribu- 
tion function with respect to conformal time, dr = dt/a, 
r is the decay rate, a is the cosmological scale fac- 
tor, qDDM is the comoving momentum, and Eddm = 
Viddm + M 2 a 2 . Specializing to two-body decays, one 
can show that the corresponding change to the SDM dis- 
tribution function will be [20T l53j 

fsDMipSDM) = 7TB ^ r / dEf DDM {p), 

4-UiSDM PSDM PCM J Ei 

(9) 

where 



Ef,i = 2 E SDMml±p S DMPcMM/m 2 SDM , 



4 



the quantity pcm is the center-of-mass momentum, and 
ml = M 2 + m 2 . 

We define the average distribution function, f®(q,r), 
and the perturbation to the distribution function, 
^i{S,q,r), for each species of particle labeled by i, ac- 
cording to 



fi(x,q,T) =f°(q,T)[l + *i(x,<l,T)} 



(10) 



where i can be the DDM, SDM, or L. Since DDM parti- 
cles are non-relativistic, their zero order phase-space dis- 
tribution is the Maxwell-Boltzmann function. The zero 
order phase-space distribution function of SDM is [20l.l25j 

fo,sDAi(q,a) = 77^^7-7; exp(-n g ) Q(ap C M-q) (11) 

where q is the comoving momentum of the SDM particle, 
a ' — q/PCMi and t q = t(a'). This can be derived from 
the fact that the decay always generates SDM particles 
with the same physical momentum pcm- In the SDM 
distribution function, the spectrum of different momenta 
arises from decays at different times, designated by the 
cosmic scale factor a' so that q = pcm®' ■ The Heavi- 
side step function Q(apcM — q) (see Eq. [11]) enforces a 
cut-off q max = apcM at a given redshift a. This maxi- 
mum momentum stems from the fact that the maximum 
momentum at a given redshift is from decay processes 
happening at that time, while SDM with lower momenta 
are from the earlier decays. To be explicit, the average 
comoving number density of SDM particles is the integral 
of fo over momentum space, 

nsDM = J q 2 dqdQ fas DM (q) ( 12 ) 

-> dn SDM (q) = / Lirq 2 dq f 0iS DM(q) (13) 
Thus fas DM can be written as 

dn S DM(q) dn S DM 1 dnsDM 



fo,SDM{q) 



fo,SDM(q) — 



q 2 dq q 2 pcMda' H(a')q 3 dt' 

(14) 

1 d(p DDM a' 3 ) 



MH(a')q 3 dt' 



This then implies Eq. (|11[) after enforcing the maximum 
momentum at g max = &PCM- 

The evolution equations for the mean energy densi- 
ties in the two dark matter components are given by the 
integrals of Eq. [8] and Eq. [9] using the unperturbed dis- 
tribution function. They read 



PDDM + 3 — pDDM — —aFpDDM 

a 



(16) 



and 



PSDM + o — (psDM + PSDM) = L n ; To PDDM (17) 



respectively. Given the DDM energy density, the decay 
product energy density pd — Psdm + P l ca n be obtained 
using the first law of thermodynamics [5J, [55j from 



da 3 pd 
dr 



-P, 



da 3 d(a 3 p DDM ) 



dr 



dr 



(18) 



This implies that the energy density evolution of the 
massless daughter particle L is 

a, , a a(M 2 -m 2 ) 

PL + 3-(p L + PL) = PL+ 4:~PL = r — — p DDM 

a a 2M Z 

(19) 



B. Perturbations 

To compute the contemporary lensing power spectra, 
it is necessary to compute the perturbations to the dark 
matter distributions and the metric. Our treatment of 
perturbations follows the conventions established in Ma 
and Bertschinger [5oj |. We will present our results in 
the synchronous gauge, because this choice lends itself 
to numerical evaluation. In the synchronous gauge, the 
Fourier transform of the Boltzmann equation can be writ- 
ten 



dln/p . h + &n ft 2 _ 1 (df 
dln7 ^--2-(*-") J 

(20) 

The DDM perturbation equations are the same as the 
well-known equations describing CDM (see Ref. |HB|), so 
we will not describe them any further. The term on the 
right-hand side of Eq. (1201) is the so-called collisional term 
representing the change in the distribution functions due 
to interactions (decays in our case). In the absence of 
non-gravitational interactions, (df/dr)c — 0. For the 
decay products, the collision terms are 



/ dfsDM_\ 
V dr J, 



2ME 



/o,ddm(1 + ^ddm) 



(21) 



and 



Iir " 



a{M 2 - m 2 )T 
2ME 



fa 



DDM 



(1 + ^ddm). (22) 



The factors m§/(2M 2 ) and (M 2 - m 2 )/(2M 2 ) that ap- 
pear in the SDM and L collision terms can be easily un- 
derstood. Consider a two-body decay in the rest frame 
of the DDM particle, A — » B + C, with corresponding 
masses rriA-, tub, and mc- The energies of B and C in 
the rest frame of A are Eg = {iti 2 a + m 2 B — m 2 c ) j (2rriA) 



and Ec = (rn 



frig) /(2m a). So these factors 



represent the ratios of energy that have been deposited 
into different daughter particle species. 

The perturbations for the massless relativistic daugh- 
ter particles may be treated in a manner analogous to 
that of massless neutrinos, save for the peculiar distribu- 
tion function of the L. We integrate out the momentum 
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dependence in the distribution function by defining (in 
Fourier space) 

P ( r ~ \ Jq 2 dqqf L (q,r)^ L {k,q,n,r) 
F L (k,n,r) = r — w 23 

An expansion of Fl in a series of Legendre polynomials 
Pi(k ■ n) has the form 

F L (k,n,r) = ^o(-i) l {'2l + l)F Ltl (k,r)P l (k-h). (24) 

The standard synchronous gauge perturbations in den- 
sity, velocity, and anisotropic stress are 



^ = ^-2^)-ar- J —-J 1 , (34) 



5*2 <? fc / 0lTf o,t, \ I 1 A i 2 ■ \ d fsDM,0 

_ = -(2*, _ 3* 3 ) - (-ft + ^)^- 

ar^^$ 2 , (35) 



fsDM,0 



and 

g*j 
9r 



(2Z 



^-—(/^-(I+l)*^!))-^ 



DDM.Q 



and 



(25) 
(26) 

(27) 



respectively. Evaluating the Boltzmann equation for our 
Legendre polynomial expansion in Eq. (1241) yields the 
evolution of the multipole coefficients in the conventional 
notation, 

Sl = + 29 L ) + a r§^5M(5 DDM - &,), (28) 

3 M pl 



9 L = k (— -at) - ar— — — 6> L , (29) 



a L = ^-(26 L + A + ffi)-H*J^)-ar§^a tl (30) 

15 4 M pL 



fsDM.O 

(36) 

for Z > 3. 

If we restrict attention only to cases in which the mass 
difference between the DDM and SDM particles is small, 
/ = 1 — m/M <C 1, the SDM particle will receive an 
extremely non-relativistic kick velocity vu — fc. As we 
should expect, SDM behaves similarly to CDM, aside 
from the fact that it is endowed with a non-negligible 
distribution of momentum due to the DDM decays. In 
this limit, the SDM perturbations evolve as for a stan- 
dard non-relativistic dark matter species, 



Ssdm 
and 

Qsdm - 
where 



-0s 



DM 



1 •• r Eip DDM 

ai — (oddm — osdm) 



-h- 



M PSDM 



(37) 



% , SPsDM 

'SDM 

SDM 



Sp 



k 2 S, 



SDM' 



r Ei pddm 

-0,1 tt-t OSDM, 



M psD 



M 



5P ; 



SDM 



and 



Fz, _ 



f a- 4 Jq 2 dq^f (q)^ 

SpsDM 47ra- 4 J q 2 dqEf {q)^ 



(38) 



(39) 
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j[lF Lil -i-{l+l)F w+1 }-aT 



Fa pddm 
M pl 



Fl,i, I > 3 
(31) 



Here we have defined E 1 = (M 2 +m 2 )/(2M) = mg/(2M) 
and E 2 = (M 2 — m 2 )/(2M). 

The SDM must be treated differently to account for 
their finite mass and non-trivial velocity kicks. We ex- 
pand the perturbation to the distribution function, *, in 
a Legendre series 

oo 

*(^n J g 1 r)=5^(-*) , (2J + l)*,(^g 1 r)fl(fc.fi). (32) 

1=0 

We have dropped the "SDM" subscript on * for brevity 
as there should be no cause for confusion in this con- 
text. Evaluating the Boltzmann evolution equation on 
this expansion, we obtain for the different multipoles 

9*o gfc l;d\D.fsDMfi. r E i Iddm ,o , t , 

"5— = E7*l + 7rft Jj har ^"? -*BCM,0 

or E 6 ding E Jsdm.q 

-ar— #o, (33) 

& JSDMfl 



The higher multipole terms become negligible in the non- 
relativistic as they are proportional to powers of the ratio 
of the kinetic energy to the total energy, q/e. 

Though we solve the complete equations for the evolu- 
tion of the SDM perturbations, the non-relativistic kick 
velocity approximation is valid in most of our calcula- 
tions. The most interesting constraints from future sur- 
veys are relevant for models with Vk ^ 10~ 3 c and rela- 
tivistic kicks have already been ruled out for a wide range 
of lifetimes (23.157} . 

Ferturbation growth is suppressed on scales smaller 
than the free- streaming scale. The free-streaming scale 
is, in turn, determined by an integral of the sound speed 
c s . We defer a detailed discussion of the free-streaming 
scale in our decaying dark matter models and its imprint 
on the matter and lensing power spectra to § IVII 



IV. NONLINEAR CORRECTIONS TO 
STRUCTURE GROWTH 

Our most robust constraints stem from perturbations 
on linear scales. However, it is interesting to estimate 
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the level of constraints that may be achieved by ex- 
ploiting mildly nonlinear scales as is common practice 
in the established framework for exploring dark energy 
with lensing and galaxy clustering statistics [3(| . Includ- 
ing mildly nonlinear scales improves constraints because 
it increases the signal-to-noise of lensing measurements 
and because it includes information regarding the effects 
of DDM on the abundance and internal structures of 
cluster-sized dark matter halos. We explore constraints 
including mildly nonlinear scales as a means of estimat- 
ing the level of constraints that may be achievable after 
an exhaustive numerical simulation program, similar to 
what is being performed for dark energy !58j . 

We implement the nonlinear corrections to the matter 
and lensing power spectra using the halo model ;59j . The 
halo model is known to exhibit mild systematic offsets 
compared to numerical simulations and the nonlinear cor- 
rection of Ref. [6(| ■ However, we use the halo model be- 
cause it provides a convenient framework for estimating 
the alterations to nonlinear structure induced by DDM 
before performing an exhaustive numerical investigation. 
We combine the standard aspects of the halo model with 
an analytical model proposed by Sanchez-Salcedo [H[ for 
the alterations to dark matter halo structure due to the 
kick velocities generated in the decay process. 

For relevant lifetimes (T -1 > Hq 1 ), dark matter halos 
begin with the same density profiles as in the standard 
CDM model. Their density distributions can be well de- 
scribed by Navarro et al. [6l[ (NFW) profiles. As the 
DDM decays, the kinetic energy of dark matter particles 
will change because SDM particles receive a small kick 
velocity from their parent particles. Assuming that we 
only consider decay processes with / « 1, the mass of 
the parent and daughter particles will be nearly iden- 
tical. As discussed in Sanchez-Salcedo jisj . on average 
the net effect of decays is to impart an amount of en- 
ergy AE m mv\j2 on the dark matter, independent of 
the initial velocity. The changes in average kinetic en- 
ergy will result in changes in particle orbits, causing an 
expansion of dark matter halos and a shallowing of dark 
matter profiles. 

To demonstrate the effect of density profile modifica- 
tion, we adopt a two-step calculation. Assume that DDM 
particles in halos follow circular orbits prior to any sig- 
nificant DDM decays. The particles orbit in the gravita- 
tional potential of the NFW halo, which can be approx- 
imately described by a power law v c (r) = vo(r /ro) 1 / 2 ^ 
over any sufficiently small range of r. In a given time in- 
terval, a small fraction of DDM particles decay and their 
daughter SDM particles gain a small amount of energy 
AE ss mvl/2. In general, the daughter particles will 
move from circular orbits to elongated orbits, character- 
ized by the new energy relative to the halo potential and 
an apocentric radius r. Orbits in the NFW potential are 
not closed, rendering it a numerical problem to compute 
the time-averaged value of the radial coordinate of the 
daughter particle. To obtain a simplistic estimate of the 
new radii the particles move to, we assume that the new 
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FIG. 1: Dark matter density profiles times radius, r p(r), as 
a function of radius and time. The dark matter halo mass 
is Mh = 1O 12 M0. In the absence of dark matter decays, 
the halo concentratio n is c = 5. The halo has a virial speed 
v-oir = sj GMh/ Rvir ~ 130 km/s. Different panels are for 
different choices of kick velocity and lifetime as labeled along 
the top and right axes respectively. In each panel the solid 
lines show the initial NFW profile. The short-dashed line, 
long-dashed line, dash-dotted line, and dash-double-doted line 
represent density profiles after 2.5, 5, 7.5, and 10 Gyr. This 
figure is designed to be directly comparable to the simulation 
results displayed in Fig. 1 of Ref. [15|] . 

average position of the daughter is similar to the radius 
of a circular orbit at the new value of the orbital energy. 
This is conservative in the present context, because cir- 
cular orbits in equilibrium are least susceptible to such 
expansion 18]. In this assumption, the radial position of 
the daughter particles, r', will be 

"-(• j/ ' + 5?Tt( ! ^)>)'- m 

The model we have described is not self-consistent, so 
it is important to validate the basic predictions of the 
model against more complete calculations. To check the 
validity of this model, we compare our analytical calcu- 
lation results with the N-body simulation results from 
[l5| . In Figure [TJ we plot density profiles for a dark mat- 
ter halo with mass Mh = 1O 12 M and an initial NFW 
concentration parameter c = 5 for several different life- 
times and kick velocities. Peter et al. computed the 
profiles of dark matter halos in the same model using 
N-body simulations that accounted for the dark matter 
decays. Fig. [T] above is the same as Figure 1 in Peter 
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FIG. 2: Similar to Figure [T] but for halos with M h = 5 x 
10 13 Mq and NFW concentration c=5. The halo virial speed 
is v V ir ~ 477 km/s. 

et al. [lH save for the fact that we have computed mod- 
ified halo profiles according to the analytic model de- 
scribed in this section. A comparison of the two figures 
reveals that the analytic model and the numerical simu- 
lations arc in remarkable agreement for all models with 
Vk ^ 200 km/s and T _1 > 10 Gyr. There are several pos- 
sible explanations for the inconsistencies that arise when 
Vk > 200 km/s and T" 1 < 10 Gyr. One is that when 
changes to the gravitational potential are not small, the 
final gravitational potential is sufficiently different from 
the initial gravitational potential that the initial poten- 
tial cannot be used to approximate the new positions of 
the SDM particles. Another possibility is that typical cir- 
cular orbits no longer provide useful approximations for 
the degree of halo expansion. As discussed in [l5[ , where 
they look at velocity anisotropy of their simulated halos, 
they found that the orbits become radially biased at the 
halo outskirts. Moreover, vj. = 200 km/s is considerable 
compared to the virial velocity of a Mh = 10 12 M Q halo, 
so it is not surprising that those halos are not in dynam- 
ical equilibrium for large Vk and small lifetime. These 
simulation results show that the assumptions of our sim- 
ple model are violated in the regime of high kick velocity 
and low lifetime. As we show in § IVI1 our primary re- 
sults in which the nonlinear model is used correspond to 
Vk < 200 km/s and lifetimes > 100 Gyr, so our use 
of this model for a first foray into this regime is justified. 

Unlike Peter et al. [15J, we are interested in cos- 
mological weak lensing as our observable. The halo 
mass most relevant to weak lensing lie in the range 



Af h w 10 13 - 10 14 M Q HI]. Such halos have significantly 
larger virial velocities than the 10 12 M halos consid- 
ered above. Typical virial velocities of these larger halos 
lie in the range v V i r s» 280 — 600 km/s. This suggests 
that our model can be used at larger Vk than the value 
Vk ~ 200 km/s that we arrived at by comparing to sim- 
ulations of a 10 12 M Q halo above, because these kicks 
represent a smaller fraction of the potential well depth. 
For instance, Peter et al. [H[ pointed out that the cluster 
mass function is insensitive to Vk ^ 500 km/s, because 
the typical virial speeds clusters are Vk ^ 600 km/s. For 
completeness, we show the corresponding density profile 
modifications for these group- and cluster-sized halos in 
Figure [U We will show in § I VII that our calculations are 
only sensitive to DDM parameters that result in density 
profiles with mild changes. 

We include this effect in our nonlinear halo model cal- 
culation by giving all recomputing halo profiles as de- 
scribed above prior to computing lensing correlations. 
We modified halo profiles by assuming initial halos with 
the same profiles, including concentrations, as their con- 
cordance ACDM counterparts and implementing the 
above model on these halos. The remainder of our halo 
model implementation follows the approach we used in 
Ref. [li], so we do not repeat it. Ideally, one would treat 
nonlinear corrections to structure growth using program 
of cosmological numerical simulations. However, we place 
such a study outside the scope of the present work as 
our initial aim is to estimate the constraining power of 
forthcoming surveys. In this manner, we estimate the 
fruit that a computationally-intensive numerical simula- 
tion program may bear on the problem of unstable dark 
matter. 



V. FORECASTING METHODS 

The Fisher Information Matrix provides a simple es- 
timate of the parameter covariance given data of spec- 
ified quality. The Fisher matrix has been utilized in 
numerous, similar contexts in the cosmology literature 
[2rj l38l lia . [62H69I , so we give only a brief review of im- 
portant results and the caveats in our particular applica- 
tion. We have confirmed the validity of the Fisher matrix 
approximation in models of unstable dark matter using 
Monte Carlo methods as described in Ref. (26| . 

The Fisher matrix of observables in Eq. (Q}, subject to 
covariance as in Eq. l[Tj). can be written as 

(41) 

where the indices A and B run over all ./Vtom^tom + 
l)/2 spectra and cross spectra, the pi are the parame- 
ters of the model, f s k y is the fraction of the sky imaged 

— 1/2 

by the experiment, and £ m in = 2/ sfcy is the smallest 
multipole constrained by the experiment. FF is a prior 
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Fisher matrix incorporating previous knowledge of viable 
regions of parameter space. We set t max = 300 for lin- 



ear forecasts and £„ 



3000 in our most ambitious 



nonlinear forecasts. On smaller scales (higher £), vari- 
ous assumptions, such as the Gaussianity of the lensing 
field, break down [38|, |48-52, 70]. To be conservative, we 
explore modest priors on each parameter independently, 
so that = 5ij/(a[) 2 , where af is the la prior on 
parameter^. The forecast, la, marginalized constraint 
on parameter pi is a{pi) = ^[F^ 1 ].^. 

Our DDM model has two independent parameters, 
namely decay rate T (or lifetime, T -1 ) and mass loss 
fraction / (which is related to Vk via Vk = fc). Either 
one of these parameters can independently be tuned to 
render the effects of DDM negligible, so it is not useful 
to marginalize over one parameter to derive constraints 
on the other. In what follows, we choose to illustrate the 
effectiveness of lensing to constrain DDM by fixing life- 
time and quoting possible constraints on /. Other than 
the mass loss fraction /, we also consider six cosmolog- 
ical parameters that we expect to modify weak lensing 
power spectra at significant levels and to exhibit partial 
degeneracy with our model parameters. We construct 
our forecasts for DDM lifetime bounds after marginaliz- 
ing over the remaining parameters. Our six additional 
parameters and their fiducial values (in parentheses) are 
the dark energy density (0.74), the present-day dark 
matter density, wdm = ^dm/i 2 (0.11), the baryon den- 
sity Ub — £lbh 2 (0.023), tilt parameter n s (0.963), the nat- 
ural logarithm of the primordial curvature perturbation 
normalization ln(A^) (—19.94), and the sum of the neu- 
trino masses J2i m ^i (0-05 eV). This choice of fiducial 
model implies a small-scale, low-redshift power spectrum 
normalization of as — 0.82. The optical depth to reion- 
ization has a negligible effect on the lensing spectra on 
scales of interest, so we do not vary it in our analysis. 

We take priors on our cosmological parameters of 
a(uj m ) = 0.007, a{uj b ) = 1.2 x 10~ 3 , a(\nA 2 R ) = 0.1, 
a{n s ) = 0.015, and a{VL&) — 0.03. We assume no pri- 
ors on DDM model parameters or neutrino mass. Our 
fiducial model is motivated by the WMAP seven-year 
result and our priors represent marginalized uncertain- 
ties on these parameters based on the WMAP seven-year 
data These priors are very conservative and allow 
for weaker constraints on DDM than would be expected 
from future data, where stronger priors may be available. 
To estimate the potential power of lensing constraints on 
DDM when stronger cosmological constraints are avail- 
able, we also explore prior constraints on these parame- 
ters at the level expected from the Planck mission 6 us- 
ing the entire Planck prior Fisher matrix of Ref. [7l| . 
Of course, using published priors from other analyses is 
not self-consistent because these priors were derived in 
analyses that assume stable dark matter, but for rele- 



vant lifetimes the dark matter decays should cause only 
subtle alterations to the cosmic microwave background 
anisotropy spectrum so this analysis should approximate 
a self-consistent simultaneous analysis of all data. 

In some cases, we will estimate nonlinear power spec- 
tra in models with significant neutrino masses. In such 
cases, we follow the empirical prescription established in 
previous studies (e.g., Refs. [33|, [3J, l68| ) and take 



where 



/„ y/P}**(k) + fb+DM ^P b N + L DM (k) 



/b+DM — 



DM 



■fife 



ill r, 



(42) 

(43a) 
(43b) 



P„{k) is the linear power spectrum of neutrinos, and 
^b+biiW i s tne nonlinear power spectrum evaluated for 
baryons and dark matter only. However, we note that 
recent work has questioned the robustness of this treat- 
ment of neutrino mass using direct numerical simulations 
[72l ] and perturbation theory 73] , so it may become nec- 
essary to revisit this aspect of the modeling of power 
spectra prior to the availability of observational data. 

We explore possible constraints from a variety of forth- 
coming data sets. We consider the Dark Energy Sur- 
vey (DES) as a near-term imaging survey that could 
provide requisite data for this test. We model DES 
by taking a fractional sky coverage of f s t~ y — 0.12 and 
with n = 15/arcmin 2 . Second, we consider a class of 
future "Wide" surveys as may be carried out by the 
Large Synoptic Survey Telescope (LSST)[27l or Euclid 
[28| . We model these Wide surveys with f s k y — 0.5 and 
n = 50/arcmin 2 . Lastly, we consider a comparably nar- 
row, deep imaging survey. We refer to such a survey 
as a "Deep" survey and model it with f s ^ y — 0.05 and 
n = 100/arcmin 2 . Such a survey may be more typical 
of a space-based mission similar to the proposed Wide- 
Field InfraRed Survey Telescope (WFIRST). In all cases, 
we take v / (t 2 T = 0- 2 and assume particular shape mea- 
surement errors from each experiment are encapsulated 
in their effective number densities, in accord with recent 
conventional practice in this regard. Our results are rel- 
atively insensitive to number density because shot noise 
does not dominate cosmic variance on the scales we con- 
sider, and our linear constraints are completely insensi- 
tive to the choice of galaxy number density over a wide 
range. 



VI. RESULTS 



1 http://www.esa.int/planck 



There are several effects of DDM on lensing power 
spectra at low redshift. First, decays change the cosmo- 
logical energy density. This change alters both structure 
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FIG. 3: Free-streaming scale as a function of scale factor. 
The blue lines show free-streaming scales for lifetime much 
greater than the age of universe ( > 100 Gyr) for several 
different mass loss fractions. The dash-dotted magenta lines 
are for / = 10~ 2 and three lifetimes. From top-to-bottom 
at right, these are 0.01 Gyr, 0.1 Gyr, and 1 Gyr. The green 
line is the free-streaming scale for massive neutrino with m^; 
=0.4 eV. Structure grows on scales between the free-streaming 
scale and horizon. On scales smaller than free-streaming scale 
(k > kps), structure growth is suppressed. 



which can be defined as 

kFs(a) = \ ~ — 7-( j 44 
V 2 c s {a) 

where %{a) = ada/dr and Ti.^ 1 is the comoving horizon 
scale. 

We show the evolution of free-streaming scale of SDM 
particles as a function of scale factor in Figure |3j for 
several mass loss fractions / and lifetimes. As dis- 
cussed in [25|, the behavior of the free-streaming scale 
of DDM can be divided into two regimes. When the 
decay process is still occurring, corresponding to cosmo- 
logical times less than the decay lifetime, daughter parti- 
cles with the same physical momentum are continuously 
created so that the sound speed stays approximately 
the same. In this case, the evolution of free-streaming 
scale will simply trace the evolution of horizon. If de- 
cays have ceased, which will happen when < Hq 1 , 
the sound speed will decrease as c s oc a -1 . The free- 
streaming scale shrinks as the initial velocities are red- 
shifted away. This effect also happens to massive neutri- 
nos as they become non-relativistic. At early times the 
neutrino free-streaming scale traces the horizon so long as 
the neutrinos have relativistic velocities. In Figure [3] we 
can see that after neutrinos become non-relativistic, at 
a nr ~ 1.3 x 10~ 3 (0.4eV/m„), their free-streaming scale 
varies as kps fll ^ 2 during matter domination, which is 
identical to free-streaming in the small lifetime limit of 
DDM. 



growth and distance. Further, decaying dark matter re- 
sults in significant free-streaming of daughter SDM par- 
ticles. While each of these effects can be important, for 
models near the limit of what may be constrained by lens- 
ing surveys, it is the effect of free-streaming that largely 
determines the lensing power spectra. In this case, 
the free-streaming velocity of SDM suppresses structure 
growth on scales smaller than free-streaming scale, an 
effect similar to that caused by massive neutrinos. In 
§ IVI Al we show the behavior of the free-streaming scale 
in our model and compare it with massive neutrinos. We 
follow this by showing the alterations to lensing power 
spectra in § IVI Bl Our final results are the forecast con- 
straints on DDM, which we give in § IVI CI 



A. Free-streaming Scale 

In the standard cosmological scenario, matter density 
fluctuations at a particular scale grow once the scale en- 
ters the horizon (k > H) during the matter-dominated 
epoch. However, species with non-negligible primordial 
velocities will be able to escape the potential wells and 
suppress the formation of structure. The scale that cor- 
responds to this effect is the free-streaming scale kps. 



B. Weak lensing Power Spectrum 

As we mentioned above, DDM affects lensing power 
spectra in two respects. First, the power spectra for po- 
tential and density fluctuations are modified by the free 
streaming of the daughter SDM particles. At k > kps, 
structure growth is suppressed. Second, the matter den- 
sity is reduced as decays occur, slightly suppressing the 
late-time growth of structure. In the left panel in Fig- 
ure[H we show that significant decrements in power occur 
at roughly the same scale, k > lCP 2 /iMpc -1 for a vari- 
ety of lifetimes, so long as the lifetime » H^ 1 (the 
regime most relevant to our work). This suppression is 
due to free streaming and indeed, the scale on which the 
suppression occurs agrees with the estimates of the free- 
streaming scale shown in Fig. [3J The right panel of Fig. 0] 
illustrates that the scale of suppression is determined by 
the mass-loss fraction /, in the limit that ^> Hq 1 . 
In models with larger /, the velocities of the daughter 
SDM particles are higher, so at fixed lifetime, they free- 
stream greater distances. Both panels in Fig. 2] show a 
small increment in power on large scales for models with 
small lifetimes (r _1 < 50 Gyr) and larger mass-loss frac- 
tions (/ > 0.1). This delineates the parameter regime for 
which the overall change in the energy budget begins to 
have a non-negligible effect on fluctuation growth. The 
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FIG. 4: Fractional difference between matter power spectrum for standard ACDM and a decaying dark matter model evaluated 
at z = 0. Left: The effect of varying the DDM lifetime at fixed mass-loss fraction, / = 1CP 1 . Solid curves show the linear 
theory predictions, and dash or dash-dot lines show predictions that include the nonlinear corrections implemented via the halo 
model. The green lines show the spectrum in a ACDM with massive neutrinos, Sm„ =0.5 eV, for comparison. Right: The 
effect of varying mass-loss fraction /, at a fixed lifetime of F _1 = 50 Gyr. 



small increment on large scales in these cases enforces a 
fixed observed CMB normalization. 

Notice in the left panel of Fig. g] that with / — 10 _1 , 
the free-streaming suppression is similar to that induced 
by massive neutrinos with the sum of the neutrino masses 
T,m u ~ 0.5 eV. This suggests that neutrinos may be de- 
generate with DDM, and this would be the case if it were 
not possible to probe a wide range of length scales and 
redshifts. In practice, we find that massive neutrinos are 
distinguishable from DDM for two reasons. First, the 
differences in scale dependence exhibited in Fig. [4] give 
a possible handle with which to separate the two. More 
importantly, the redshift dependence of the power spec- 
trum differs in the two models. This is most easily seen 
in Fig. [3j The evolution of the free-streaming scale of 
massive neutrinos and the free-streaming scale of DDM 
differs significantly. Deep, large-scale survey data that 
enable probes of structure at a variety of redshifts be- 
tween < z < 3, as is expected of forthcoming surveys, 
break the potential degeneracy between massive neutri- 
nos and DDM. 

The observed strength of gravitational lensing also has 
a dependence upon geometry, so differences in angular di- 
ameter distance may lead to modified lensing power spec- 
tra. Changes in relative partitioning of energy among 
relativistic and non-relativistic species will change the 
evolution of the angular diameter distance. However, 
we consider small mass loss fractions (/ <Cl) and large 
lifetimes (T _1 >> H^ 1 ), so angular diameter distances 
are altered only by negligible amounts and, although we 
account for these changes, they do not provide lever- 
age on constraining DDM. For example, for a lifetime 
of T -1 = 50 Gyr and / = 10 _1 , Figure g] shows that 



the changes in matter power spectrum are approximately 
40-60% for k > 2 x 10~ 2 /iMpc -1 ). The corresponding 
changes in angular diameter distances are less than 0.1%. 
The constraining power of weak lensing comes primarily 
from the suppression of structure growth. Incidentally, 
this is a promising feature because the primary informa- 
tion used to constrain dark energy using lensing surveys 
is carried by the geometric piece of the lensing signal 

MM- 

Aoyama et al. [25j have considered constraints on un- 
stable dark matter stemming from contemporary data on 
the cosmological distance-redshift relations. Using con- 
straints on the Hubble parameter, the baryon acoustic 
oscillation scale, and the angular positions of the cos- 
mic microwave background anisotropy spectrum peaks, 
they have placed competitive constraints on such mod- 
els. Aoyama et al. [25| nn d that T" 1 > 0.1 Gyr with 
/ < 3.5 x 10" 2 and that T" 1 > 30 Gyr when / ~ 1. 

The Dash-dotted lines in Figure [5] exemplify the al- 
terations to the small-scale lensing convergence power 
spectra incurred when we account for the altered halo 
profiles that result from dark matter decays. As / in- 
creases, kick velocities increase, and the fractional power 
decrement increases, as we should expect. This addi- 
tional suppression is confined to relatively small scales 
(large multipoles, £ > 300) for most of the parameter 
space of interest (v k < 200 km/s for T" 1 < 100 Gyr). 

As we pointed out in right panel in Figure HI DDM 
may partially mimic massive neutrinos if redshift evolu- 
tion information in not accessible. In Figure |6l we show a 
comparison of the redshift evolution of DDM and massive 
neutrino lensing power spectra in three tomographic red- 
shift bins. Other than the difference in shapes, it is also 
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FIG. 5: Fractional difference of auto convergence lensing 
power spectrum between standard ACDM model and decay- 
ing dark matter model from first tomographic redshift bin 
(lensing source galaxies between < z v < 0.6, where z v is 
photometric, and not necessarily true redshift). Solid lines 
are calculated using halo model with NFW profiles. These 
lines include the alteration of the linear power spectrum on 
large-scales and the reduction in the abundance of dark mat- 
ter halos due to free-streaming. However, halos are assumed 
to have the same profiles as they would in standard ACDM. 
The Dash-dotted lines include the nonlinear corrections to 
halo density profiles. 



evident that the DDM power spectra evolve significantly 
more than the spectra in massive neutrino models. The 
reason is that the decay process continuously deposits ki- 
netic energy into the daughter dark matter distribution, 
in contrast to the neutrinos which have purely redshifting 
kinetic energy distributions. 



C. Forecast Constraints on DDM Model 
Parameters 

To estimate of the power of weak lensing to constrain 
DDM, we adopt a variety of possible strategies. First, 
we consider constraints from data on scales where linear 
evolution of density fluctuations should be valid. The 
value of this approach is that exploiting linear scales to 
constrain DDM does not require a simulation program 
to confirm or refine nonlinear models of structure for- 
mation in these models. This can be done with con- 
temporary theoretical knowledge of the phenomenology 
of these models. Moreover, relatively large-scale con- 
straints are less observationally challenging because they 
exploit data on scales where cosmic variance, rather than 
galaxy shape measurements, are the dominant error [2fj| . 
In both cases, these constraints are conservative so we 
should expect that forthcoming lensing surveys designed 



FIG. 6: Comparison of the redshift evolution of decaying 
dark matter and massive neutrino lensing power spectra. We 
plot fractional difference of auto convergence lensing power 
spectra between standard ACDM model and decaying dark 
matter (or massive neutrino) models in three tomographic 
redshift bins (labeled at the top) . For simplicity, we show only 
the linear power spectra in this plot, though spectra computed 
with our nonlinear model lead to a similar conclusion. 

to address dark energy should do at least as well as our 
linear forecasts. To limit ourselves to linear scales, we 
take data on multipoles £ < 300. All of the constraints 
that we show in this section have been marginalized over 
the remaining cosmological parameters, including neu- 
trino mass. 

To show the maximum potential of lensing surveys, 
we consider measurements that extend into the mildly 
nonlinear regime, as is commonly done for dark energy 
forecasts. The primary value of this extension is not that 
particular features in the power spectra induced by DDM 
are added to the data set. Rather the primary improve- 
ment in constraints comes from an increase in the signal- 
to-noise with which the power suppression can be de- 
tected 26]. In this case, we include information on multi- 
poles up to our quoted maximum multipole £ max = 3000 
(see § [Vj. Constraints on these scales will rely on reli- 
able modeling of clustering on mildly nonlinear scales, 
so a comprehensive simulation program will be necessary 
to ensure the robustness of such constraints. A com- 
prehensive program is computationally-intensive and be- 
yond the scope of our present paper, as part of our goal is 
to emphasize that such a large-scale numerical program 
may be interesting and useful. 

In Figure [7] we display our forecast la exclusion con- 
tours alongside a variety of other contemporary con- 
straints. The most relevant contemporary constraints 
come from modifications to the structures of dark matter 
halos with virial velocities similar to the SDM kick veloc- 
ities [l5| (orange region) . Additional constraints may be 
placed on unstable dark matter by examining the proper- 
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FIG. 7: Comparison of DDM parameter exclusion contours from Peter et al. flEI ] (orange) and Peter and Benson [53] (dark 
and light green) to those that may be achieved from weak gravitational lensing. The exclusions in |.22j rely on a variety of 
assumptions regarding the merger history and formation histories of galaxies. Roughly speaking these parameter regimes also 
correspond to parameter values for which the interpretation of the missing satellites problem may be significantly altered by 
unstable dark matter (see Ref. f22| for details). The red, blue, and purple lines are the 1 a exclusion contours from our weak 
lensing forecasts assuming "wide," "deep," and DES like galaxy imaging surveys respectively. The solid lines combine weak 
lensing with projected Planck constraints, and the dot-dash lines are from weak lensing alone. The less restrictive set of contour 
lines correspond to weak lensing constraints using scales on which linear theory is applicable (£ < 300). The more restrictive 
set of contours incorporate multipoles up to i" m ax = 3000 and the theoretical calculation includes our nonlinear corrections to 
halo density profiles. 



ties of the satellite galaxies of the Milky Way (22| (green 
regions). However, these constraints rely on a variety 
of assumptions regarding the formation and evolution of 
relatively small galaxies. Moreover, these constraints de- 
lineate a range of DDM parameters for which unstable 
dark matter may have a significant effect on the inter- 
pretation of the missing satellites problem. As this type 
of DDM model has been invoked to mitigate the " miss- 
ing satellite" problem, it should not be a surprise that 
stronger constraints also come from these types of obser- 
vations. As such, it is this parameter range for which 
it is most interesting to develop independent constraints 
on unstable dark matter and that is the purpose of our 
weak lensing study. 

As indicated in Figure our most conservative, lin- 
ear calculation can already give interesting constraints 
DDM that are competitive with contemporary bounds. 
The largest advantage of lensing constraints will be that 
it can extend constraints on DDM lifetimes significantly, 
as is evident in Fig. [7J These forecasts are not depen- 
dent upon modeling nonlinear structure growth, so they 



constitute a robust lower limit to the constraining power 
of imaging surveys. Moreover, these constraints are not 
subject to any particular assumptions regarding the evo- 
lution and formation of galaxies, particularly the Milky 
Way satellite galaxies that are the subject of so much 
contemporary research. Comparing the linear constraints 
from the three types of surveys, the "Deep" survey pro- 
vides slightly more restrictive constraints than DES. A 
"Wide" survey similar to LSST or Euclid has the poten- 
tial to improve the constraints relative to DES by ~40- 
60%. 

The slope of the constraint contours turns over near 
lifetimes of a few Gyr. This turn over reflects the turn 
over in free streaming scale exhibited by the dash-dotted 
lines in Fig. [3] In models with -C Hq , the free- 
streaming scale at low rcdshift decreases with time. No- 
tice that including Planck priors yields only a marginal 
improvement on the forecast constraints, ~ 15 — 40% over 
the parameter ranges of interest. Our nonlinear forecasts 
exhibit a similar sensitivity to Planck priors, though they 
are not depicted in Fig.[7Jin the interest of clarity. 
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Our more ambitious approach is to estimate con- 
straints including data on mildly nonlinear scales. As 
we have already mentioned, forecast constraints in this 
case are less robust because they rely on data in a regime 
on which systematic errors and galaxy shape noise, as op- 
posed to cosmic variance, will likely be the largest errors 
and because theoretical models of structure growth on 
these scales will need to be validated using a large suite of 
cosmological structure formation simulations. The most 
ambitious constraint forecasts in Fig. [7] result from in- 
cluding multipoles I < £ max — 3000 along with Planck 
priors on the standard cosmological parameters. These 
forecasts indicate that weak lensing constraints on Vk 
that include mildly nonlinear scales may improve upon 
the minimal, linear constraints by nearly an order of mag- 
nitude. In this case, it is clear that weak lensing may pro- 
vide an independent and competitive constraint on DDM 
and that lensing constraints will extend to significantly 
longer lifetimes than contemporary bounds. 



VII. DISCUSSION AND CONCLUSION 

In this paper we have explored the effect of decaying 
dark matter on the large-scale matter distribution and 
the power of future weak lensing surveys to place con- 
straints on DDM lifetimes and mass splittings. The mass 
difference is parameterized by the velocity kick (vk) that 
the daughter, stable dark matter particles receive upon 
the decay of the heavier, parent DDM. DDM leads to 
a suppression of matter clustering on scales below the 
free-streaming scale of the stable, daughter dark matter 
particles and this suppression can be probed with data 
from galaxy imaging surveys. 

Our most conservative constraint forecasts result from 
considering lensing over large scales on which linear the- 
ory should be valid. In this case, best limits which come 
from a "Wide" survey, similar to Euclid or LSST. These 
surveys may exclude Vk > 90 km/s for ^1 — 5 Gyr, a 
result that is competitive with contemporary constraints 
[l5j |. Lensing improves upon contemporary constraints 
most markedly for large decay lifetimes (r^ 1 > H^ 1 ). In 
this regime, lensing constraints are significantly more re- 
strictive than contemporary bounds from halo structure. 
In the relatively near-term, the DES will be able to place 
limits of Vk ^ 160 km/s for ^1 — 5 Gyr. Achieving 
constraints at this level should be achievable. First, the 
lensing surveys we study are under development to study 
dark energy already. Moreover, these constraints assume 
that we restrict attention only to relatively large scales 
on which linear perturbation theory can be used to pre- 
dict lensing power, so no additional theoretical effort will 
be necessary. 

It may be possible to derive more restrictive lensing 
constraints on unstable dark matter by considering the 
mildly nonlinear scales that are commonly considered as 
part of the program to constrain dark energy. Including 
multipoles up to ell ~ 10 3 increases constraining power 



by boosting the signal-to-noise with of the weak lens- 
ing signal on scales that are sensitive to the dynamics of 
the dark matter. Exploiting such scales will rely on an 
exhaustive simulation program to understand nonlinear 
clustering in DDM models, similar to the simulation pro- 
gram that is being performed in support of dark energy 
probes 58] , so significant additional theoretical work will 
be necessary. Nevertheless, the payoff may be significant. 
In order to estimate the ambitious constraints that may 
be achieved from such a data analysis, we have imple- 
mented nonlinear corrections to lensing power using the 
standard halo model coupled with a simple model for the 
modification of halo density structures due to decaying 
dark matter. 



In our most ambitious forecasts, we find that weak 
lensing may constrain the mass splitting of the DDM 
nearly an order of magnitude more restrictively than im- 
plied by our linear scale analysis. These results sug- 
gest that weak lensing surveys will be sensitive to Vk ~ 
10 km/s for lifetimes T _1 < 10 Gyr for all the survey 
types that we have considered. These constraints are in- 
teresting because they restrict parameters for which con- 
siderable effects due to DDM may be evident in the Milky 
Way satellite galaxy population. It may be possible to 
achieve similar constraints depending upon a variety of 
assumptions regarding the formation process of these 
satellite galaxies (22| , but lensing provides a complemen- 
tary constraint using data on distinct length scales. 



We have demonstrated that measurements of the large- 
scale matter distribution through a weak lensing survey 
will be a powerful probe of decaying dark matter. This 
probe is valuable for several reasons. First, such surveys 
as PanSTARRS, LSST, DES, Euclid, and WFIRST are 
already being undertaken as part of the effort to con- 
strain dark energy. The survey requirements specified 
by the dark energy program are the same that we as- 
sume here, so no additional observational work will be 
necessary. Moreover, we have shown that such measure- 
ments can provide independent, competitive constraints 
on models of DDM that could alter our interpretation 
of the small-scale problems of the standard cosmologi- 
cal model, particularly the missing satellites problem. In 
fact, we have demonstrated that lensing will probe DDM 
models with lifetimes that exceed contemporary bounds 
by an order of magnitude. Our most ambitious constraint 
forecasts rely upon the development of accurate and pre- 
cise models of matter clustering in models of unstable 
dark matter. This will likely require a significant simu- 
lation effort to ensure the robustness of any constraints 
derived from forthcoming data. It is our hope that this 
proof-of-concept work will motivate more detailed numer- 
ical studies of unstable dark matter models as well as ad- 
ditional possible constraints from related observations. 
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